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ABSTRACT 

Motivated by the recent detection by the Wilkinson Microwave Anisotropy Probe of a large optical 
depth to Thomson scattering - implying a very early reionization epoch - we assess a scenario where 
the universe was reionized by 'miniquasars' powered by intermediate- mass black holes (IMBHs), the 
remnants of the first generation of massive stars. Pregalactic IMBHs form within minihalos above the 
cosmological Jeans mass collapsing at z > 20, get incorporated through mergers into larger and larger 
systems, sink to the center owing to dynamical friction, and accrete cold material. The merger history 
of dark halos and associated IMBHs is followed by Monte Carlo realizations of the merger hierarchy in 
a ACDM cosmology. Our model is based on the assumptions that quasar activity is driven by major 
mergers and nuclear IMBHs accrete at the Eddington rate a fraction of the gas in the merger remnant. 
The long dynamical frictional timescales leave many IMBHs 'wandering' in galaxy halos after a minor 
merger. While seed IMBHs that are as rare as the 3.5cr peaks of the primordial density field evolve 
largely in isolation, a significant number of black hole binary systems will form if IMBHs populate the 
more numerous 3cr peaks instead. In the case of rapid binary coalescence a fraction of IMBHs will 
be displaced from galaxy centers and ejected into the IGM by the 'gravitational rocket' effect, rather 
than accrete and shine as miniquasars. We show that, under a number of plausible assumptions for the 
amount of gas accreted onto IMBHs and their emission spectrum, miniquasars powered by IMBHs - and 
not their stellar progenitors - may be responsible for cosmological reionization at z ~ 15. Reionization 
by miniquasars with a hard spectrum may be more 'economical' than stellar reionization, as soft X- 
rays escape more easily from the dense sites of star formation and travel further than EUV radiation. 
Energetic photons will permeate the universe more uniformly, make the low-density diffuse IGM warm 
and weakly ionized prior to the epoch of reionization breakthrough, set an entropy floor, and reduce gas 
clumping. Future 21 cm observations may detect a preheated, weakly-ionized IGM in emission against 
the cosmic microwave background. 
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1. INTRODUCTION 



In popular cold dark matter (CDM) hierarchical cos- 
mogonies, the very first generation of metal-free ('Popula- 
tion III') stars is expected to form in dark matter 'miniha- 
los' of total mass M h > 5 x 10 5 M (e.g. Fuller & Couch- 
man 2000; Yoshida et al. 2003) condensing due to H2 cool- 
ing from the high-cr peaks of the primordial density field 
at redshift z — 20 — 30. Recent numerical simulations 
of the collapse and fragmentation of primordial molecu- 
lar clouds suggest that the first stars were predominantly 
very massive, m* > 100 M Q , a scale linked to H 2 chem- 
istry and cooling (Bromm, Coppi, & Larson 2001; Abel, 
Bryan, & Norman 2002). Zero-metallicity very massive 
stars have spectra similar to that of a ~ 10 5 K blackbody 
(Tumlinson & Shull 2000), and emit about 20 times more 
Lyman-continuum photons per stellar baryon than a stan- 
dard stellar population (Schaerer 2002; Bromm, Kudritzki, 
& Loeb 2001). A 'top-heavy' primordial initial mass func- 
tion (IMF), quite different from the present-day Galactic 
case, may have then played a crucial role in determining 



the ionization, thermal, and chemical enrichment history 
of the intergalactic medium (IGM) at early times. 

The IGM is known to be highly ionized at least out to 
redshift ~ 5.5. While the excess H 1 absorption measured 
in the spectra of z ~ 6 quasars in the Sloan Digital Sky 
Survey (SDSS) has been interpreted as the signature of the 
trailing edge of the cosmic reionization epoch (Djorgovski 
et al. 2001; Fan et al. 2002; White et al. 2003), the recent 
detection by the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite of a large optical depth to Thomson 
scattering, r e = 0.17 ± 0.04 (68%), suggests that the uni- 
verse was reionized at much higher redshift, z ro j = 17 ± 3 
(assuming instantaneous hydrogen reionization, see Kogut 
et al. 2003; Spergel et al. 2003). This is an indication of 
significant star- formation activity at very early times. In 
the wake of WMAP remarkable results, several theoretical 
studies of stellar reionization scenarios have pointed out 
that an early generation of massive, metal-free stars with 
a top-heavy IMF may be necessary to produce a Thom- 
son optical depth as high as the central value determined 
by WMAP (Cen 2003a; Ciardi, Ferrara, & White 2003; 
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Haiman & Holder 2003; Sokasian et al. 2003; Wyithe & 
Loeb 2003; Somerville & Livio 2003). 

Since, at zero metallicity, mass loss through radiatively- 
driven stellar winds or nuclear-powered stellar pulsations is 
expected to be negligible (Kudritzki 2002; Baraffe, Heger, 
& Woosley 2001), Population III stars will likely die losing 
only a small fraction of their mass. Nonrotating very mas- 
sive stars in the mass window 140 < to* < 260 M© will 
disappear as pair-instability supernovae (Bond, Arnctt, & 
Carr 1984), leaving no compact remnants and polluting 
the universe with the first heavy elements (e.g. Schnei- 
der et al. 2002; Oh et al. 2001; Wasserburg & Qian 2000). 
Stars with 40 < to* < 140 M and to* > 260 M are pre- 
dicted instead to collapse to black holes (BHs) with masses 
exceeding half of the initial stellar mass (Heger & Woosley 
2002). Barring any fine tuning of the IMF of Population 
III stars, intermediate-mass black holes (IMBHs) - with 
masses above the 4-18 M Q range of known 'stellar-mass' 
BHs (e.g. McClintock & Remillard 2003) - may then be 
the inevitable endproduct of the first episodes of pregalac- 
tic star formation. Since they form in high-cr density fluc- 
tuations, 'relic' IMBHs and their descendants will tend to 
cluster in the bulges of present-day galaxies as their host 
halos aggregate into more massive systems (Madau & Rees 
2001, hereafter Paper I; Islam, Taylor, & Silk 2003). The 
first IMBHs may also have seeded the hierarchical assem- 
bly of the supermassive variety (SMBHs) observed at the 
center of luminous galaxies (Volontcri, Haardt, & Madau 
2003, hereafter Paper II; Volonteri, Madau, & Haardt 
2003, hereafter Paper III) . 

Population III IMBHs accreting gas from the surround- 
ing medium will shine as 'miniquasars' at z <~ 15 and 
generate a soft X-ray background that may catalyze the 
formation of H2 molecules in dense regions and counter- 
act the destruction by UV Lyman- Werner radiation (Cen 
2003b; Glover & Brandt 2003; Haiman, Abel, & Rees 
2000). The net effect would be an increase in the cool- 
ing rate and star formation efficiency of minihalos (but 
see Machacek, Bryan, & Abel 2003 who find the positive 
feedback effect of X-rays to be quite mild). In this pa- 
per we point out that miniquasars represent an additional 
source of Lyman-continuum photons that must be taken 
into account in models of early reionization by Population 
III objects. Thin disk accretion onto a Schwarzschild black 
hole releases about 100 MeV per baryon. If 'seed' IMBHs 
were able to (say) double their initial mass via gas accre- 
tion, and just a few percent of the radiated energy were 
emitted above and close to the hydrogen Lyman edge, then 
miniquasars would be more efficient at photoionizing the 
universe than their metal-free stellar progenitors. 

While the basic idea is simply stated, it is not easy to 
quantify and model in detail. A minihalo at z ~ 20 will 
typically undergo several major mergers in a Hubble time. 
Pregalactic IMBHs will follow the merger history of their 
hosts, and may be able to accrete efficiently only in the in- 
ner densest regions of the merger remnant. Owing to the 
long dynamical frictional timescales, many IMBHs will not 
sink to the center and will be left 'wandering' in galaxy 
halos after a minor merger. If seed holes are numerous 
enough, binary systems may form in significant numbers. 
In this case, rather than accrete and shine as miniquasars, 
a fraction of IMBHs may be displaced from galaxy cen- 
ters and ejected into the IGM by the 'gravitational rocket' 



effect (Redmount & Rees 1989) or by triple BH interac- 
tions (Paper II). Miniquasars are expected to be copious 
sources of soft X-ray photons, which will permeate the 
IGM more uniformly than possible with extreme ultravio- 
let (EUV, > 13.6 eV) radiation (Oh 2000, 2001) and make 
it warm and weakly ionized prior to the epoch of reioniza- 
tion breakthrough (Venkatesan, Giroux, & Shull 2001). 

In this paper we make a first attempt at treating the 
growth and impact on the very early IGM of Population III 
IMBHs in the context of hierarchical structure formation 
theories, incorporating some of the essential astrophysics 
into a scenario for the reheating and reionization of the 
universe by miniquasars. The plan is as follows. In § 2 we 
outline a model in which seed holes populate the rare 3.5cr 
peaks of the primordial density field. We show that, if cold 
material can be accreted efficiently onto IMBHs hosted in 
minihalos, then miniquasars may be responsible for cos- 
mological reionization at z <~ 15. In § 3 we present results 
for a case in which seed holes are more numerous but gas 
accretion occurs less efficiently, and address the dynam- 
ics of IMBH binary systems in shallow potential wells. In 
§ 4 we discuss the implications of an early reionization 
and reheating epoch by hard radiation from miniquasars. 
Finally, we summarize our results in § 5. Unless other- 
wise stated, all results shown below refer to the currently 
favoured ACDM world model with Q, M = 0.3, fi A = 0.7, 
h = 0.7, n b = 0.045, cr 8 = 0.93, and n = 1. 

2. BUILD UP OF IMBHs IN MERGING MINIHALOS 

A scenario for the birth, growth, and evolution of 
IMBHs that traces their formation in minihalos with virial 
temperatures r v ; r < 10 4 K back to very high redshifts was 
put forward in Papers I and II. Seed IMBHs are placed 
far up in the dark halo merger tree, get incorporated into 
larger and larger halos, sink to the center owing to dynami- 
cal friction, and accrete a fraction of the gas in the merger 
remnant. The merger history of DM halos and associ- 
ated IMBHs is followed through Monte Carlo realizations 
(based on the extended Press-Schechter formalism) of the 
merger hierarchy from early times to the present. 

Here we improve on several aspects and adapt some pre- 
scriptions in our model to the conditions expected dur- 
ing the epoch of first light. Pregalactic seed IMBHs 
form within the mass ranges 20 < to. < 70 M Q and 
130 < to. < 200 M©, as remnants of the first generation 
of massive metal-free stars that do not disappear as pair- 
instability supernovae. For simplicity, within these inter- 
vals the differential black hole mass function is assumed to 
be flat, i.e. the IMF-averaged BH mass is 115 M Q . They 
form in isolation within minihalos above the cosmological 
Jeans mass collapsing at z = 24 from rare > v-a peaks of 
the primordial density field. As our fiducial model we take 
v = 3.5, corresponding in the assumed ACDM cosmology 
to minihalos of mass M seod = 1.3 x 10 6 M Q . As halos more 
massive than the v-a peaks contain a fraction erfc(f/\/2) 
(= 0.00047 for v = 3.5) of the mass of the universe, the 
mass density parameter of our '3.5<r' pregalactic holes is 

n. = °- 00047 ^ (w - ) = 10-™ (i) 

M sccd 

This is 0.6% of the mass density of the supermassive vari- 
ety found in the nuclei of most nearby galaxies, Osmbh = 
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(2.1 ± 0.3) x 1(T 6 (Yu & Tremaine 2002). We also 
study a case in which seed holes are more numerous and 
populate all peaks above 3a. Then erfc(3/v / 2)=0.0027, 
Msccd = 1-2 x 10 5 Mq, and the initial mass density of '3<t' 
IMBHs is fi. = 10 -6 - 1 , less than 40% of the mass density 
of SMBHs today. We do not consider lighter, more com- 
mon halos than those originating from 3a fluctuations, 
as they would fall below the minimum mass threshold 
for baryonic condensation via molecular hydrogen cool- 
ing, A/th = 1-25 x 10 5 Mq according to the simulations 
of Machacek, Bryan, & Abel (2001). Seed holes that were 
much rarer than the 3.5a density peaks would have little 
impact on the reionization of the IGM at z > 15. Also, the 
assumed 'bias' assures that almost all halos above 10 11 M Q 
actually host a BH at later epochs. For simplicity, we start 
our calculations by placing all seed IMBHs in their host 
halos at z — 24, the highest redshift we are able to follow 
the merger hierarchy to because of computational costs. 
More realistically, IMBHs are expected to form continu- 
ously from Population III massive stars over a range of 
redshifts. 

Each halo is treated as a singular isothermal sphere 
(SIS) with density p oc r~ 2 , truncated at the virial ra- 
dius. During the merger of two halos, the 'satellite' pro- 
genitor with initial mass M s - and its IMBH if it hosts 
one - spirals in toward the center of the more massive pre- 
existing system (with mass Mh > M s ) on the dynamical 
friction (against the DM background) timescale. This de- 
pends on the orbital parameters of the satellite (van den 
Bosch et al. 1999; Ghigna et al. 1998) and on tidal mass 
loss/evaporation. N-body simulations show that massive 
satellites with mass ratios P = M s /Mh > 0.1 ('major' 
mergers) sink rapidly without significant mass loss, while 
the lightest satellites (P < 0.01) are almost unaffected by 
dynamical friction and do not suffer orbital decay (Taffoni 
et al. 2003). In our calculations we adopt the expression for 
the dynamical friction timescale of a mass-losing satellite 
in a main halo suggested by Colpi, Mayer, & Governato 
(1999). 
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Fig. 1. — Fraction of BH host halos that contains two or more 
seed IMBHs as a function of redshift, averaged over 50 'trees'. Two 
curves are shown, one assuming BHs form in isolation within mini- 



halos collapsing at z = 24 from 3.5<r peaks, the other from 3cr peaks 
instead. For illustrative purposes, dynamical interactions between 
BHs have been switched off along cosmic history. 

In our fiducial model (y = 3.5), the density of seed holes 
is small enough that the merging of two minihalos both 
hosting a BH is a rare event, and IMBHs evolve largely 
in isolation. This is shown in Figure 1, where the number 
of halos containing two or more primordial IMBHs, nor- 
malized to the total number of BH hosts, is plotted as a 
function of redshift. The ratio, as expected, increases with 
time and tends to saturate at lower redshifts. For illustra- 
tive purposes, the realizations in this figure have been run 
neglecting any dynamical interaction between BHs along 
cosmic history. At the very high redshifts of interest here 
this is a good assumption for v = 3.5 as the fraction of 
multiple systems is typically less than 10%. If seed holes 
populate the v — 3 peaks instead, the fraction of multiple 
BH hosts can exceed 50% and bound binaries form in sig- 
nificant numbers. We will discuss the dynamics of IMBH 
binaries in § 3. 

2.1. Accretion onto IMBHs 

Hydrodynamic simulations of major mergers have shown 
that a significant fraction of the gas in interacting galaxies 
falls to the center of the merged system (Mihos & Hern- 
quist 1996): the cold gas may be eventually driven into 
the very inner regions, fueling an accretion episode and 
the growth of nuclear BHs (Kauffmann & Hachnclt 2000; 
Cavaliere & Vittorini 2000). Since the local SMBH mass 
density is consistent with the integrated luminosity den- 
sity of quasars (Yu & Tremaine 2002; Fabian 2003), the 
fraction of cold gas ending up in the hole must depend 
on the properties of the host halo in such a way to ul- 
timately lead to the observed correlation between stellar 
velocity dispersion and SMBH mass (Ferrarese & Merritt 
2000; Gebhardt et al. 2000). In Paper II it was assumed 
that in each major merger the hole in the most massive 
halo accretes at the Eddington rate a gas mass that scales 
with the fifth power of the circular velocity of the host 
halo: the normalization was fixed a posteriori in order to 
reproduce the observed local itibh — cr* relation (Ferrarese 
2002). Our simple model was shown to explain remarkably 
well the observed luminosity function of optically-selected 
quasars in the redshift range 1 < z < 5. Accretion was 
inhibited onto all BHs hosted in halos with virial tem- 
perature T v ; r < 10 4 K, as below this temperature atomic 
cooling is not effective in allowing the gas to sink to the 
center. 

Prior to the epoch of reionization, however, cold ma- 
terial may be efficiently accreted onto IMBHs hosted in 
minihalos. As mentioned in § 1, gas condensation in the 
first baryonic objects is possible through the formation of 
H2 molecules, which cool efficiently via roto-vibrational 
transitions even at virial temperatures of a few hundred 
kelvins. In the absence of a UV photodissociating flux 
and of ionizing X-ray radiation, three-dimensional simula- 
tions of early structure formation show that the fraction of 
cold, dense gas available for accretion onto IMBHs or star 
formation exceeds 20% above Mh = 10 6 M Q (Machacek, 
Bryan, & Abel 2003). And while radiative feedback (pho- 
todissociation and photoionization) from the progenitor 
massive star may initially quench BH accretion within the 
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original host minihalo, new cold material will be readily 
available through the hierarchical merging of small gaseous 
subunits. In this section we will assume that gas accretion 
onto BHs can actually occur with high efficiency in mini- 
halos just above the cosmological Jeans mass. 

As there is no compelling reason to expect the tobh — <?* 
relation to be set in primordial structures and be satisfied 
at the very high redshifts [z ~ 15) considered in this work, 
we run two different sets of realizations. In the first we 
assume that in each major merger the hole in the main 
halo triples its mass via gas accretion, A77z acc — c 1ttl^yv 1 
while in the second the hole accretes a fixed fraction of 
the host halo mass, Am acc = 10~ 3 Mfr,. In both real- 
izations this mass is added to the hole at the Eddington 
rate, rriE = AirGm p mBti/ '(core) , for a radiative efficiency 
e = 10%. In the case of a merger in which only the satel- 
lite halo contains a BH, the hole is placed at the center of 
the merger remnant after one dynamical friction timescale, 
and must await a subsequent major merger event with 
a (smaller) satellite to start accreting. Each accretion 
episode begins after one dynamical timescale (measured 
at radius 0.1r v ; r ) and typically lasts for a few Salpeter 
timescales, t$ = wbh/wb = 4.5 x 10 7 yr (e/O.l) -1 , until 
a mass Am acc has been added to the hole. Figure 2 shows 
how the comoving mass density of IMBHs builds up with 
cosmic time: with the adopted prescriptions, by a redshift 
of 14 this never exceeds 3% of the mass density of SMBHs 
inferred by Yu & Tremaine (2002) in local galaxies. 
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Fig. 2. — Comoving mass density in IMBHs as a function of red- 
shift for the two different accretion recipes discussed in the text. Gas 
accretion onto BHs is assumed to occur in minihalos only after a ma- 
jor merger. Seed IMBHs form in minihalos collapsing at z = 24 from 
3.5<t density peaks. 



2.2. Miniquasars 

Assume now that a fraction /uv of the bolometric power 
radiated by our Population III miniquasars is emitted as 
hydrogen- ionizing photons with mean energy (hv). For a 
given bolometric emissivity, the number of ionizing pho- 
tons scales as /uv/C 1 ^)- One of the biggest uncertain- 
ties in discussing early reionization by miniquasars is their 
unknown emission spectrum. If the shape of the emit- 



ted spectrum from miniquasars followed the mean spec- 
tral energy distribution of the quasar sample in Elvis 
et al. (1994), then /uv = 0.3 and (hv) — 3ryd, hence 
fvv/(hv) = 0.1 ryd -1 . Miniquasars powered by IMBHs, 
however, are likely to be harder emitters than quasars. 
The hottest blackbody temperature in a Keplerian disk 
damping material onto a BH at the Eddington rate is 
fcT'max ~ lkeV?7J B ]!/ 4 (Shakura & Sunyaev 1973), where 
the hole mass is measured in solar mass units. The char- 
acteristic multicolor (cold) disk spectrum follows a power- 
law with L u ~ v 1 / 3 at hv < kT max . Close to the 
epoch of reionization at z ~ 15, most IMBHs holes will 
have masses in the range 200 — 1000 M Q (Fig. 3), hence 
ff max ~ 0.2 — 0.3 keV. The spectra of 'ultraluminous' X- 
ray sources (ULXs, Colbert & Mushotzky 1999) in nearby 
galaxies actually require both a soft and a hard compo- 
nent of comparable luminosities to describe the continuum 
emission. While the soft components are well fit by cool 
multicolor disk blackbodies with kT max ~ 0.15 keV, which 
may indicate IMBHs (e.g. Miller et al. 2003), the nonther- 
mal power-law component has spectral slope L v cx v~ a , 
with a « 1. A power-law with slope a — 1 extending from 
2 keV down to 13.6 cV has /uv = 1 and (hv) = 5ryd, 
hence fw/(hv) = 0.2 ryd . Figure 4 shows the cumu- 
lative number of ionizing photons produced per hydrogen 
atom for different values of /uv / (hv) . A few ionizing hard 
photons per atom (less than for a stellar spectrum because 
of secondary ionizations and larger mean free path, see § 4) 
should suffice to reionize the universe and keep the gas in 
overdense regions and filaments photoionized against ra- 
diative recombinations. 
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Fig. 3. — Mass function of accreting IMBHs predicted at four dif- 
ferent redshifts by our fiducial model. Units are arbitrary. All seed 
holes are assumed to form at z = 24 from 3.5cr density peaks. 
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Fig. 4. — Cumulative number of ionizing photons per hydrogen 
atom produced by miniquasars, for different values of the fraction 
/uv of the bolomeric power that is emitted as radiation above 1 ryd 
with mean energy (hu) (in ryd). Solid curves: in each major merger 
the BH in the main halo accretes a mass Am aC c = 2rrtBH- Dashed 
curves: same with Am aC c = 10 -3 M h . Dotted curve: same with 
Am acc = 10 — 3 M^, fuv/ (hi/) = 0.2 ryd -1 , but with gas accretion 
suppressed in minihalos with virial temperature T v i r < 10 4 K. 

To illustrate the implications of these results, consider 
the following estimate for the number of H-ionizing pho- 
tons emitted by the initial population of progenitor mas- 
sive stars. In our model the fraction, /*, of cosmic baryons 
incorporated into Population III massive stars - progeni- 
tors of seed IMBHs - at z = 24 is, 

/* = ^^5xl0- (2) 

where fib is the baryon density parameter. Zero- 
metallicity stars in the range 40 < m* < 500 M© emit 
about 70,000 photons above 1 ryd per stellar baryon 
(Schaerer 2002). The total number of ionizing photons 
produced per baryon by progenitor Population III stars is 
then, using equation (1), 

Tl' 

w 70, 000/» w 0.04, (3) 

n b 

well below what is needed to reionize the universe. Figure 
4 shows that, if fuv /(his) is greater than 0.1 ryd -1 and 
gas is accreted efficiently onto IMBHs, then miniquasars 
may be responsible for cosmological reionization at redshift 
~ 15. We will discuss in more details some of the implica- 
tions of reheating by hard radiation from miniquasars in 
§4. 

3. REDUCED ACCRETION EFFICIENCY AND DYNAMICS 
OF IMBH BINARIES 

Miniquasars powered by IMBHs forming in 3.5<7 peaks 
produce more ionizing photons than their progenitors, and 
will reionize the IGM by z ~ 15. Such conclusion is 
based on the assumption that IMBHs can grow rapidly 
via gas accretion during major mergers involving miniha- 
los and shine in the EUV/soft X-rays. Local and global 



feedback effects, however, may act to reduce or even in- 
hibit accretion in minihalos where atomic cooling is in- 
efficient. More massive halos with T V j r > 10 4 K, in which 
atomic cooling can operate, are expected to be very rare at 
these redshifts. If gas accretion was completely suppressed 
in minihalos with T v ; r < 10 4 K, i.e. with masses below 
10 7 - 9 [(l+z)/20]~ 3 / 2 M Q , the lower emissivity would either 
shift the reionization epoch to lower redshift or make mini- 
quasars a sub-dominant source of ionizing photons (Fig. 
4), at least until halos massive enough that gas can ini- 
tially cool and contract via excitation of hydrogen Lya, 
collapsed in large numbers. 

Still, a viable scenario for early reionization by mini- 
quasars may be one in which gas accretion onto IMBHs is 
much reduced compared to our fiducial model in § 2, but 
seed holes are more numerous at the start in order to sus- 
tain the early production of ionizing radiation. Black holes 
forming in 3cr density fluctuations cannot be assumed to 
evolve in isolation, and the dynamics of BH binaries in 
shallow potential wells must be addressed. In this section 
we shall assess a model in which IMBHs populate all halos 
more massive than M see( j = 1.2 x 10 5 M©, correspond- 
ing to the 3<T peaks of the primordial density field. In 
each major merger a mass Am acc = 10~ 5 Mh (two orders 
of magnitude smaller than in our fiducial model) is now 
added to the hole in the more massive halo at the Edding- 
ton rate, i.e. seed IMBHs do not grow appreciably due 
to gas accretion and miniquasars are typically 'on' only 
for a small fraction of the Salpeter timescale. Note that 
a high accretion recipe with Am acc = 10~ 3 Mh (say) is 
not viable in the case of 3c seed holes, as it generates a 
Pbh that exceeds the mass density of SMBHs observed in 
nearby galaxies. 



3.1. IMBH binaries 

As shown in Figure 5, if seed holes are as numerous as 
the 3cr density peaks, binary systems may form in signifi- 
cant numbers. The merging - driven by dynamical friction 
against the DM - of two minihalo+IMBH systems with 
mass ratio P > 0.1 will drag in the smaller hole towards 
the center, leading to the formation of a bound IMBH 
binary in the violently relaxed core of the newly merged 
system. In massive galaxies at low redshift, the subse- 
quent evolution of the binary may be largely determined 
by the central stellar distribution. The binary initially 
shrink ('hardens') by dynamical friction from distant stars 
acting on each BH individually, and later via three-body 
interactions, i.e., by capturing the stars that pass within 
a distance of the order of the binary semi-major axis and 
ejecting them at much higher velocities (Begelman, Bland- 
ford, & Rees 1980). Dark matter particles will be ejected 
by decaying binaries in the same way as the stars, i.e. 
through the gravitational slingshot. Their contribution to 
binary hardening will be approximately weighted by the 
fraction /dm = Pbm/(pdm + p*), typically <C 1 in massive 
galaxy cores. 
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Fig. 5. — Fraction of all halos hosting nuclear black holes (NBHs) 
that contain a hard BH binary (BHB) system, as a function of red- 
shift. Binaries are 'hard' if their separation is a < (see text for 
details). Seed IMBHs form in minihalos collapsing at z = 24 from 3<j 
density peaks. In this figure we have assumed that black hole bina- 
ries do not shrink and coalesce (i.e. we set the hardening timescale 
t h -» oo). 



In minihalos a numerous population of low-mass stars 
may be present if the IMF were bimodal, with a second 
peak at 1-2 M , as suggested by Nakamura & Umemura 
(2001). Otherwise the binary will be losing orbital energy 
to the dark matter background. In the limit /dm 3> 1 (cf. 
Paper III), a binary with BH masses mi > ni2 and semi- 
major axis a(t) in an SIS dark halo of density p and one- 
dimensional velocity dispersion a, becomes 'hard' when 
the binary separation falls below the value 



ah 



Gm 2 
4cr 2 



0.04 pc 



m 2 



10 3 M 



(4) 



(Quinlan 1996), where 05 is measured in units of 5 kms . 
The hardening of the binary modifies the density profile, 
removing mass interior to the binary orbit, depleting the 
galaxy core of dark matter, and slowing down further hard- 
ening. If .M C j is the mass ejected by the BH pair, the 
binary evolution and its effect on galaxy cores are deter- 
mined by two dimensionless quantities: the hardening rate 



H 



a da 1 
Gp dt 



and the mass ejection rate 



J 



dM ci 



(mi + 1JI2) d\na 1 



(5) 



(6) 



The quantities H and J can be found from scattering ex- 
periments that treat the test particle-binary encounters 
one at a time (Quinlan 1996). We assume as in Papers II 
and III that the removal of matter creates a core of radius 
r c and constant density p c = p(r c ), eroding a preexisting 



isothermal cusp with p oc 



The total mass ejected as 



the binary shrinks from ah to a is then M e j = 4cr 2 r c /3G, 
while the core radius grows as 



3 f ah 
r c(t) = —~G(mi + m 2 ) / J(a)da/a 

4cr Ja(t) 



(7) 



The binary separation quickly falls below r c and subse- 
quent evolution is slowed down due to the declining matter 
density, with a hardening time, 



t h = \a/a\ 



2irr c (t) 2 
Haa 



(8) 



that becomes increasingly long as the binary shrinks. If 
hardening continues down to a separation 



1 AU 



(mi 



m-2 )mim2 



10 9 - 3 M„ 3 



1/4 



(9) 



the binary will coalesce within 200 Myr (the Hubble time 
at z = 18) due to the emission of gravitational waves. The 
pair must then manage to shrink by a factor a,h/a gv ~ 5000 
for gravity wave emission to become efficient. Since the 
hardening and mass ejection rate coefficient are H as 15 
and J ps 1 in the limit of a very hard binary (Quinlan 
1996), the hardening time in equation (8) can be rewrit- 
ten (for mi = m 2 ) as 



t h w 1.3 x 10 5 yr 



m 2 



10 3 M m 



af^- ln 2 (a h /a). (10) 



This becomes longer than the then Hubble time already 
at ah I a ~ 100. Note that for an equal mass binary 

aft/ag r cx m^a^ 2 . Thus even if the tobh oc <j 4 relation 
(Tremaine et al. 2002) were to hold at very high redshifts, 
for gravitational radiation to induce coalescence a hard 
BH binary must decay by a larger factor in minihalos with 
small velocity dispersion than in more massive galaxies. 
At these early epochs, when the BH occupation fraction 
is still small, we find that triple BH interactions arc too 
rare to cause binaries to shrink (cf. Paper II). In the ab- 
sence of any other mechanism for removing orbital angular 
momentum, the decay is then expected to stall at a sep- 
aration much greater than a gr . This is even more true as 
it is the total matter density that was allowed to decrease 
in equation (7) , not the density of low-angular momentum 
particles that get close enough to extract energy from the 
binary, i.e. we have assumed the 'loss cone' is constantly 
refilled. Note that this would be a good assumption if or- 
bital angular momentum losses to stars were the dominant 
mechanism of binary shrinkage, since the stellar relaxation 
timescale 



0.34cr 3 



G 2 m*p c In A 



150 Myr 



0.5 pc 



In A 



^5, 



is typically much shorter than the hardening time in mini- 
halos, and two-body scatterings will efficiently refill the 
loss cone. 

3.2. Gravitational rocket and the depletion of IMBHs 
from host minihalos 

The extent with which BH pairs loose angular momen- 
tum to DM particles or stars is one of the major uncertain- 
ties in computing merger timescales, and makes it difficult 
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to construct viable merger scenarios for BH binaries (e.g. 
Milosavljevic & Merritt 2001). For illustrative purposes, 
we run here two different sets of realizations. In the first 
the hardening is due to the ejection of DM particles, the 
hardening timescale is long (from eqs. 8-10), and bina- 
ries stall. Another possibility is that gas processes, rather 
than three-body interactions with stars or DM, may in- 
duce IMBH binaries to shrink rapidly and coalesce (e.g. 
Armitage & Natarajan 2002; Gould & Rix 2000). Gas in a 
system containing a hard binary will fall inward in just the 
same way as if there were a single hole at the centre, as the 
binary separation is much smaller than the Bondi radius, 
r B = Gm BH A4 = 0.5pc(m BH /1000 M Q ) (T/1000K)" 1 . 
It will heat up and radiate away its binding energy. Some 
material will be accreted by each individual hole while 
some will be flung out, removing binary orbital angu- 
lar momentum and hastening the merger process. Fur- 
ther infall may replenish the gas reservoir on a dynamical 
timescale and keep the drag going. The details are un- 
clear. As an extreme case we have therefore run a sec- 
ond set of realizations in which once a hard binary forms 
with separation a/,, it coalesces instantaneously. In the 
shallow potential wells of minihalos, the growth of IMBHs 
will then be halted by the 'gravitational rocket', the recoil 
due to the non-zero net linear momentum carried away 
by gravitational waves in the coalescence of two unequal 
mass black holes. Radiation recoil is a strong field ef- 
fect that depends on the lack of symmetry in the sys- 
tem, and may eject IMBHs from the cores of minihalos. 
To date, the outcome of a gravitational rocket remains 
uncertain, as fully general relativistic numerical compu- 
tations of radiation reaction effects are not available at 
the moment. The predicted recoil velocity, vqm, may be 
bracketed (Brandt & Anninos 1999; Wiseman 1992) on 
the one side by quasi-Newtonian calculations predicting 
vcm < 100 kms" 1 (Fitchett 1983), and on the other by 
extrapolated perturbative results on test particles (Fitch- 
ett & Detweiler 1984), allowing for rocket velocities up to 
vcm — 700 kms" 1 for mi = 1.2m2 in a Schwarzschild 
geometry. 

According to quasi-Newtonian calculations for circular 
orbits, while the binary shrinks its center of mass spirals 
outward with a velocity 



«cm = 1500 kms 



f(m 1 /m 2 ) 



(12) 



where r s = 2G(mi + ni2)/c 2 is the Schwarzschild radius 
of the system, r ms = 6Gtoi/c 2 is the radius of the closest 
stable circular orbit for two non-rotating holes (Clark & 
Eardley 1977), and the function 



/(m 1 /m 2 ) = |(l + ^i)- 2 (l + ^)- 3 (l 
m-2 mi 



m 2 

771l 



)l, (13) 



can be the most disruptive and lead to the larger modifica- 
tions in terms of IMBH occupation fraction. The coalesced 
binary will leave the galaxy altogether if its recoil velocity 
exceeds the escape speed of the halo, giving origin to a 
population of intergalactic IMBHs. 6 




14 16 18 20 22 1 
redshift 



4 16 18 20 22 24 
redshift 



Fig. 6. — Contribution of nuclear, wandering, and intergalactic 
holes to the IMBH mass density as a function of redshift. In each 
major merger a mass Am acc = 10 — 5 is accreted onto the BH 
in the main halo. Left panel: model in which (most) binaries stall. 
Right panel: model in which binaries shrink rapidly. Solid line: 
mass density of IMBHs in galaxy nuclei. Long-dashed line: wander- 
ing IMBHs retained in galaxy halos, mostly due to minor mergers. 
Short-dashed line: intergalactic IMBHs ejected from the host after 
a gravitational rocket. The horizontal dotted line shows the mass 
density of SMBHs in the nuclei of nearby galaxies inferred by Yu & 
Tremaine (2002). 



Figure 6 shows the predicted nuclear, wandering, and 
intergalactic IMBH mass density as a function of redshift, 
for the two sets of realizations described above. Two fea- 
tures are worth noting, the actual decrease with time of 
the density of nuclear IMBHs as more and more holes be- 
come 'wandering' after minor mergers (cf. Fig. 2), and 
the one order of magnitude increase in the density of in- 
tergalactic holes due to the gravitational rocket after rapid 
coalescence. In the case where (most) binaries stall, the 
small population (in mass density) of intergalactic IMBHs 
is due to a few mergers involving near-equal mass light 
binaries in the most massive halos, where the hardening 
timescale (eq. 10) is shorter. In the model where binaries 
coalesce rapidly we find recoil velocities that exceed the 
escape speed from their hosts in 80% of the cases. Indeed, 
because of the rocket effect, the mass density of intergalac- 
tic holes exceeds that in the nuclear variety below a redshift 
of 15 or so. 



reaches a maximum value / max = 0.01789 for mi = 2.6777,2. 
The extrapolated perturbative results (Fitchett & De- 
tweiler 1984) yield considerably higher recoil velocities for 
777,1 < 4 7772, while for m\ > 1 mi both calculations con- 
verge to recoil velocities below 15 kms" 1 . We adopt here 
the extrapolated perturbative results, as this assumption 

6 We have assumed that the final mass ttibh of the hole after coalescence satisfies the entropy-area relation (maximally efficient radiative 
merging): nijj = + m|. 
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Fig. 7. — Cumulative number of ionizing photons per hydrogen 
atom produced by miniquasars for different values of /tjv / (hu) (in 
units of ryd -1 ), assuming Am acc = 10 — 5 M^. Seed IMBHs form in 
minihalos collapsing at z = 24 from 3<r density peaks. Dot-dashed 
curves: model in which BH binaries stall. Solid curves: same with 
rapid binary coalescence and gravitational rocket. 

The ejection into the IGM of IMBHs by the gravitational 
rocket effect appears - depending on the initial density of 
seed holes and the ability of BH binaries to shrink and 
coalesce - to decrease only slightly (by about 15%) the 
'reionization efficiency' of miniquasars. This can be read- 
ily seen in Figure 7, which shows the cumulative num- 
ber of H-ionizing photons produced by miniquasars for 
the two cases ('binary stalling' vs. 'instantaneous coa- 
lescence') discussed in this section. Even with a large re- 
duction of the mass accreted by each IMBH, the increased 
number of seeds produces more than a few ionizing pho- 
tons per hydrogen atom by redshift 15 for a spectrum with 
fxsv/{hv) > 0.1 ryd" 1 . If seed holes are as numerous as the 
3cr density peaks, the total number of Lyman-continuum 
photons produced per baryon by their stellar progenitors 
is now 

71 • 

— » 70,000/, « 2.5, (14) 
n b 

i.e. sixty times larger than in the 3.5a case (cf. eq. 3). 
In this case miniquasars may still dominate the ionizing 
photon budget for fuv/{hv) <; 0.2ryd _1 (or, of course, if 
Am acc > W~ 5 M h ). 

4. REIONIZATION BY HARD PHOTONS 

We have seen above that, if cold gas in minihalos can be 
accreted either efficiently onto rare ('3.5cr') Population III 
IMBHs or much less efficiently onto more numerous ('3cr') 
seed holes, then miniquasars may be responsible for cos- 
mological reionization at early times. In this section we 
discuss the reheating and reionization history of an IGM 
photoionizcd primarily by a power-law metagalactic flux 
extending to X-ray energies, rather than by photons near 
the Lyman limit, as would be the case for purely stellar 
UV radiation. Reionization by energetic photons has been 

7 The exact form of the photoionization cross-section for H I was 
(1996) for He I (see also Miralda-Escude 2000). 



discussed recently by Venkatesan, Giroux, & Shull (2001) 
and Oh (2001). 

4.1. Soft X-ray background radiation 

Accretion onto IMBHs may be an attractive way to 
reionize the low-density IGM. A large fraction of the UV 
radiation from massive stars may not escape the dense 
sites of star formation, or may be deposited locally in halo 
gas that recombines almost immediately. The harder ra- 
diation emitted from miniquasars is instead more likely to 
escape from the hosts into intergalactic space, and may 
then produce more 'durable' (albeit partial) ionization in 
the diffuse IGM. 

In the pre-reionization universe, when the IGM is pre- 
dominantly neutral, photons with energies above 24.6 eV 
will be absorbed as they photoionize hydrogen or he- 
lium atoms. For a mixture of H and He with cosmic 
abundances, the effective bound-free absorption cross- 
section 7 can be approximated (to an accuracy of 30% 
in the range 50 eV < hu < 2keV) as o^f ~ 1.6 x 
10 _17 cm 2 (hu/1 ryd) -3 . The continuum optical depth of 
a uniform IGM of hydrogen density no(l + z) 3 , for photons 
emitted at frequency u e by a source at redshift z e beyond 
the redshift of reionization z Te u is given by 



T = 



dz— n (l + z) 3 a bi {u) 



' 10 



,6.4 



'lryd^ 



(15) 



where v = u e (l+z)/(l+z e ), dl/dz ~ H^ 1 fi M 1/2 (l+^)- 5 / 2 , 
and we have evaluated the integral from z re \ — 14 to 
z e = 19. The universe is optically thick even to hard 
radiation, and all photons below (Zi^)thick ~ 2 keV will be 
absorbed across a Hubble volume. Since a spectrum with 
u L v const (like the nonthermal component observed 
in ULXs) is characterized by equal power per logarithmic 
frequency interval, photoelectric absorption by the host 
galaxies will not significantly attenuate the ionizing en- 
ergy flux. 

The mean free path of ionizing photons in the neutral 
IGM is 



A = 



1 



n (l + z) 3 cr b f 



13 pc 



l + z 
20 



hu 
1 ryd 



(16) 



At z = 24, the comoving space density of halos collapsing 
from 3.5er density peaks is ~ 15 Mpc~ 3 , corresponding to a 
mean proper distance between neighboring halos of 15 kpc. 
Photons with energies above {hu) oveT \ &v ~ 150 eV will then 
have a mean free path greater than the mean separation 
between sources. We can therefore separate the radiation 
field into two components: a fluctuating EUV field in the 
range 13.6 < hu < 150 eV, Jtjv, which creates expanding 
patchy H II regions, and a nearly uniform soft X-ray back- 
ground component in the range 150 eV < hu < 2keV, Jx, 
which partially ionizes the IGM in a fairly homogeneous 
manner. 

4.2. Secondary ionizations 

Prior to the epoch of reionization breakthrough, before 
the overlapping of fully ionized (by EUV photons) bubbles 

taken from Spitzer (1978), and we have used the fit from Verner et al. 
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around individual miniquasars, soft X-ray radiation will 
produce a warm (T IGM ~ few x 100 — 1000 K), weakly ion- 
ized (x ~ 0.1) IGM (Venkatesan et al. 2001). X-rays alone 
do not produce a fully ionized medium, but can partially 
photoionize the gas by repeated secondary ionizations. A 
primary nonthermal photoelectron of energy E = 1 keV 
in a medium with residual ionization (from the recombi- 
nation epoch) x — 2 x 10~ 4 will create over two dozens 
secondary electrons, depositing a fraction /; on w 37% of 
its initial energy as secondary ionizations of hydrogen, and 
only /heat ~ 13% as heat (Shull & Van Steenberg 1985). 




IGM ionized fraction 

Fig. 8. — Number of H secondary ionizations per ionizing photon 
vs. IGM ionized fraction. A power-law radiation field with energy 
spectrum Jx oc v~ a , extending from 150 eV to 2 keV, has been as- 
sumed. The three curves from top to bottom are for a = 0.5, 1, 1.5, 
respectively. As the medium becomes more ionized an increasing 
fraction of energy is deposited as heat. 

Figure 8 shows the mean number of secondary ioniza- 
tions per ionizing photon as a function of IGM ionized 
fraction, for a power-law spectrum extending from 150 eV 
to 2 keV ('uniform soft X-ray background'). The calcula- 
tion uses the analytic fits given by Shull & Van Steenberg 
(1985). The timescale for electron-electron encounters re- 
sulting in a fractional energy loss / = AE/E, 

(where E is measured in keV), is typically much shorter 
that the electron Compton cooling timescale off cos- 
mic microwave background (CMB) photons, tc — (7 x 
10 6 yr) [(1 + z)/20]~ 4 , and thus the primary photoelectron 
will ionize and heat the surrounding medium before it is 
cooled by the CMB. 

The history of reheating and reionization by a hard v^ 1 
spectrum may then be described as follows. Initially, when 
the IGM is mostly neutral, each photoionization produces 
a host of secondary collisional ionizations, with about one 
hydrogen secondary ionization for every 37 eV of energy 
in the primary photoelectron. For photon energies above 
0.2 keV, the ratio of the H i /He i photoionization cross- 
sections drops below 4%: since the primordial ratio of he- 



lium to hydrogen is about 8%, the photoionization and 
heating rates will be dominated by helium absorption, 
which exceeds the hydrogen contribution by ~ 2 : 1 . While 
He i is the main source of hot primary photoelectrons, 
however, H i undergoes the bulk of secondary ionizations. 
Secondary ionizations operate efficiently for H i and He i , 
but not for He n (Shull & Van Steenberg 1985). Because 
of losses to secondary ionizations and to collisional excita- 
tions of H i and He i before thermalization, the heating by 
energetic photons is initially not as efficient as for photons 
near the photoelectric thresholds. The amount of heat de- 
posited into a unit proper volume of the IGM by the mean 
field can be written in first approximation as 

AQ = /heat — / dvJ x = (1 + z) 3 f hcat f x ep acc c 2 , 

C J ^overlap 

(18) 

where fx is the fraction of the bolometric emissivity of 
miniquasars that is radiated between Overlap and i^thick, 
fx = In Khick/ ^overlap) /hi {v 2 /ui) for a v' 1 spectrum ex- 
tending from v\ to v-i- For the spectrum we have chosen, 
the result is only logarithmically sensitive to the integra- 
tion bounds, and we set fx — 0.4. In the scenarios for the 
formation and growth of IMBHs discussed in the previous 
sections, gas accretion onto IMBHs rapidly build up (by 
z ~ 22) a comoving mass density p &cc in excess of a few 
hundred M Mpc -3 , enough to heat up the gas to initial 
temperatures: 

T = J^L « 1500K ( ( ^ ) . (19) 

3/2nfc \ 0.13 J V200 M Mpc" 3 / v ' 

This same radiation field would ionize the IGM to: x w 
4 x 10~ 3 (/ion/0.37) (pacc/200 M Mpc" 3 ). Since both the 
hydrogen recombination and the Compton cooling time are 
much longer than the then Hubble time for x < 0.1, the 
gas evolves adiabatically. Once the IGM ionized fraction 
increases to x « 0.1, the number of secondary ionizations 
per ionizing photon drops to a few (Fig. 8), and the bulk 
of the primary's energy goes into heat (/heat ~ 0.6) via 
elastic Coulomb collisions with thermal electrons. At this 
point Compton cooling will restrict any further increase 
in temperature: it is the rate of direct photoionizations by 
UV photons in the soft end of the spectrum that drives 
the final stages of the reionization process, until the fully 
ionized bubbles associated with individual sources finally 
overlap. 

As H i and He i reionization nears completion, the only 
means for a hard radiation field to inject energy into 
the plasma is via He n photoionization and Compton 
heating. Compton heating is much less efficient than 
He ii photoionization heating: the small Thompson cross- 
section implies that the universe is always optically thin 
to Compton heating of electrons. Furthermore, the en- 
ergy exchange per interaction is only a small fraction of 
the photon energy, ~ hv(hv/m e c 2 ). Thus, in most sit- 
uations, Compton heating will play a sub-dominant role 
compared to He n photoheating (cf. Madau & Efstathiou 
1999). 

4.3. Entropy floor 

As recently pointed out by Oh & Haiman (2003), soft 
X-rays can exert another important effect on the IGM and 
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subsequent star formation, as the warm, partially ionized 
medium produced by energetic photons has a high entropy. 
Gas at the mean density that is heated to temperatures 

T - m>9 °K30Wk)(^o)" /3 (20) 

(where 5 is the overdensity of the gas in the minihalo in 
the absence of preheating) will have entropy in excess of 
that acquired by gravitational shock heating alone. It will 
therefore have excess pressure after adiabatic heating and 
compression. As shown in the previous section, such tem- 
peratures are easily achievable with miniquasars. This en- 
tropy floor: (a) greatly reduces gas clumping, curtailing 
the number of photons needed to maintain reionization; 
and (b) results in significantly lower gas densities in the 
cores of minihalos that suppress rapid H 2 formation. The 
latter effect may imply that X-rays inhibit rather than en- 
hance star formation (Oh & Haiman 2003). 

An entropy floor could conceivably also reduce the ef- 
ficiency with which miniquasars can shine. In halos that 
have had their central densities reduced and are unable to 
form stars, IMBHs may still accrete gas, but at a consid- 
erably reduced rate. Miniquasars, however, are ~ 15 — 150 
times more efficient at producing ionizing radiation per 
processed baryon than metal-free stars, and thus need a 
much smaller reservoir of cold gas to reionize the universe. 
If fuv/{hv) ~ 0.2, and 100 MeV are released per baryon, 
then only ~ 7 x 10~ 7 of all baryons need to be accreted 
onto IMBHs to provide one ionizing photon per baryon, 
compared to ~ 10 -5 for massive metal-free stars (Bromm, 
Kudritzki, & Loeb 2001), and ~ 10~ 4 for metal-free stars 
with a Salpeter IMF (Tumlinson & Shull 2000). The cold 
gas that cooled prior to the establishment of the entropy 
floor (and formed an accretion disk around the central 
black hole) will still be available as fuel. If only a frac- 
tion /cod ~ 1% of the gas in a 3.5ct peak halo cooled to 
form an accretion disk around the central black hole be- 
fore an entropy floor was established, this still represents 
~ 5x 10~ 6 (./cooi/0.01) of all baryons, enough to provide 
10 ionizing photons per hydrogen atom; this number in- 
creases by an order of magnitude if 3cr peaks are involved 
instead. By contrast, these baryon fractions are insuffi- 
cient to reionize the universe with Population III metal- 
free stars; the entropy floor would prohibit full reioniza- 
tion by such objects. The combination of X-rays (which 
reduce gas clumping, and thus the number of UV photons 
require to complete reionization) and UV photons (which 
carry reionization to completion) could allow miniquasars 
to reionize the universe quite efficiently. 

5. SUMMARY 

This paper should be read as a first attempt at treat- 
ing the impact on the very early IGM of Population III 
IMBHs - the remnants of the first generation of massive 
stars - in the context of hierarchical structure formation 
theories. We have tried to incorporate some of the es- 
sential astrophysics into a scenario for the reheating and 
reionization of the universe by miniquasars. In our model 
quasar activity is driven by major mergers and IMBHs 
are able to accrete at the Eddington rate only in the dens- 
est inner regions of the merger remnant. If seed IMBHs 
are as rare as the 3.5ct peaks of the primordial density 



field, they evolve largely in isolation; a significant number 
of BH binary systems forms if IMBHs populate the more 
numerous 3cr peaks instead. In the case of rapid binary 
coalescence, rather than accrete and shine as miniquasars, 
a fraction of IMBHs will be displaced from galaxy cen- 
ters and ejected into the IGM by the 'gravitational rocket' 
effect. Note that, for gravitational radiation to induce co- 
alescence, a hard BH binary must decay by a larger factor 
in minihalos with small velocity dispersion than in more 
massive galaxies. The loss of orbital angular momentum 
to a gaseous disk, rather than via three-body interactions 
with DM particles or stars, may drive IMBH binaries to 
merge rapidly, but the details are unclear. 

Under a number of different assumptions for the amount 
of gas accreted onto IMBHs and their emission spectrum, 
miniquasars powered by IMBHs - and not their metal-free 
stellar progenitors - may be responsible for cosmological 
reionization at z ~ 15. The WMAP data on r e may then 
be setting constraints on gas accretion onto the earliest 
generation of IMBHs, as well as on the IMF of the first 
burst of star formation in the universe. Even if the 'UV 
radiation efficiency' of miniquasars was lower and Popula- 
tion III stars were the dominant source of ionizing radia- 
tion at z ~ 25, the transition from Population III to Pop- 
ulation II at a redshift possibly determined by the level 
of metal enrichment of the IGM would decrease by more 
than an order of magnitude the emission rate of Lyman- 
continuum photons (e.g. Cen 2003b). Hard radiation from 
miniquasars could then sustain ionization in the redshift 
range z ro ; < z < 20 (say), before ordinary stars with a 
Salpeter IMF took over. 

At z > 4 — 5, the known population of quasars appears 
not to contribute significantly to the hydrogen ionizing 
background (Madau, Haardt, & Rees 1999). The comov- 
ing space density of luminous SDSS quasars at z ~ 6 is 
20 times smaller than that at z <~ 3 (Fan et al. 2003). 
The dearth of bright QSOs at z <~ 6 may be reconciled 
with the existence of a substantial population of accret- 
ing IMBHs at z <~ 15 if there were two distinct periods of 
quasar activity. The first may have ended with the demise 
of miniquasars due to radiative feedback mechanisms af- 
ter the reionization epoch (as photoionization from back- 
ground UV photons unbinds the gas in the shallow po- 
tential wells of minihalos, Barkana & Loeb 1999; Haiman, 
Abel, & Madau 2001) or to a raised entropy floor. Ef- 
ficient gas accretion onto nuclear BHs could conceivably 
start again only at z < 10 < z ro i, in T V1I 3> 10 K halos 
further down in the merger hierarchy that are not sus- 
ceptible to Jeans smoothing and gas cooling suppression 
effects. This may lead to a second "quasar era" seeded by 
more massive (~ 5000 M Q ) pre-existing holes. 

One might worry that with a sufficiently hard ionizing 
spectrum H i and He n— >Hc m reionization may be simul- 
taneous, whereas observations of the He n Lya forest have 
been interpreted as evidence that the double reionization 
of helium occurred at a redshift of ~ 3 (e.g. Kriss et 
al. 2001; Reimers et al. 1997). For a v~ x miniquasar 
spectrum there are 3 times more He n ionizing photons 
[hu > 4ryd) per He atom than photons above 1 ryd per 
H atom. As doubly ionized helium rccombines about 6 
times faster than hydrogen, the number of > 4 ryd pho- 
tons per He atom emitted in one recombination timescale 
is only half of the corresponding value for hydrogen, and 
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there will be just a small delay between the complete over- 
lapping of H ii and He in regions (Miralda-Escude & Rees 
1993; Madau & Meiksin 1994). He n reionization by mini- 
quasars at z > 10 would provide a significant boost to 
the gas temperature, with interesting implications for the 
thermal history of the IGM (e.g. Hui & Haiman 2003). 
The patchy He n absorption observed at z ~ 3 could then 
be explained by variations in the spectra of the ionizing 
sources (Smette et al. 2002), or may imply that He n was 
reionized twice: a first time by miniquasars, followed by 
recombination as Population II stars with softer spectra 
took over, and then by a second overlap phase at z <~ 3 
during the second quasar epoch. 

Absorption in the intervening IGM will make mini- 
quasars at z ~ 15 inaccessible to direct observations from 
2 /im down to soft X-ray energies. He n recombination 
lines longward of H i Lya may be potentially detectable by 
the James Webb Space Telescope only from the brightest 
miniquasars (Oh, Haiman, & Rees 2001). An alternative 
way to probe the end of the cosmic 'dark ages' and discrim- 
inate between different reionization histories is through 21 
cm tomography of neutral hydrogen (Madau, Meiksin, & 
Rees 1997). In general, 21 cm spectral features will dis- 
play angular structure as well as structure in redshift space 
due to inhomogeneities in the gas density field, hydrogen 
ionized fraction, and spin temperature. As recently dis- 
cussed by Tozzi et al. (2000), Ciardi & Madau (2003), and 
Furlanctto, Sokasian, & Hcrnquist (2003), the fluctuations 
in 21 cm line induced by the 'cosmic web' that develops at 
early times in CDM-dominated cosmologies could be de- 
tected in emission against the CMB provided that: (a) at 
epochs when the IGM is still mainly neutral, the first gen- 
eration of (Population III) stars produce enough UV con- 
tinuum photons with energies between 10.2 and 13.6 eV to 
mix the hyperfine levels and decouple the hydrogen spin 
temperature from Tcmb; and (b) the IGM is 'warm', i.e. 
Tigm > Tcmb = 57 K [(1 + z)/20]. X-ray radiation natu- 



rally produces a warm, weakly ionized IGM. Miniquasars 
turning on at early stages may then make the IGM visible 
in 21 cm emission as structure develops, before the uni- 
verse is actually reionized. Brightness temperature fluctu- 
ations may be detectable by future facilities like the LOw 
Frequency ARray (LOFAR) if foreground contamination 
from unresolved extragalactic point sources (Di Matteo et 
al. 2002; Oh & Mack 2003) can be successfully removed 
by using spectral structure in frequency space. 

Finally, we remark that accreting IMBHs may release 
large amounts of energy to their environment in other 
forms than radiative output. Outflows from miniquasars 
may inject kinetic energy into intergalactic space and raise 
the temperature of the IGM to a much higher adiabat than 
expected from photoionization. In the type of scenario dis- 
cussed in § 3, gas accretion along cosmic history builds up 
by z ~ 14 a comoving mass density, p acc , in excess of 
10 4 M Q Mpc -3 . If a fraction f w of the accreted rest- mass 
energy was used to drive an outflow and ultimately de- 
posited into the IGM, the energy input per baryon at this 
epoch, 

£ „ = / „^„, keV (A_) ( 104M ^ pc _ 3 ), 

(21) 

would suffice to preheat vast regions of the universe to tem- 
peratures above a few xl0 5 K, and so inhibit the forma- 
tion of early dwarf galaxies (Benson & Madau 2003). The 
detailed thermal history of the IGM may then depend on 
different forms of 'feedback' mechanisms from miniquasars 
during the epoch of galaxy formation. 

Support for this work was provided by NASA grant 
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